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Abstract 

We consider a two-dimensional model of inflation, where the inflationary trajectory is "deformed" by a 
grazing encounter with an Extra Species/Symmetry Point (ESP) after the observable cosmological scales 
have left the Hubble radius. The encounter entails a sudden production of particles, whose backreaction 
causes a bending of the trajectory and a temporary decrease in speed, both of which are sensitive to initial 
conditions. This "modulated" effect leads to an additional contribution to the curvature perturbation, which 
can be dominant if the encounter is close. We compute associated non-Gaussianities, the bispectrum and 
its scale dependence as well as the trispectrum, which are potentially detectable in many cases. In addition, 
we consider a direct modulation of the coupling to the light held at the ESP via a modulaton field, a mixed 
scenario whereby the modulaton is identified with a second inflaton, and an extended Extra Species Locus 
(ESL); all of these scenarios lead to similar additional contributions to observables. We conclude that inflaton 
interactions throughout inflation are strongly constrained if primordial non-Gaussianities remain unobserved 
in current experiments such as PLANCK. If they are observed, an ESP encounter leaves additional signatures 
on smaller scales which may be used to identify the model. 
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I. INTRODUCTION 



In the simplest models of inflation, the inflaton smoothly rolls down its potential while con- 
tinuously producing super-Hubble perturbations via amplification of quantum fluctuations. When 
trying to embed the inflaton into a more concrete framework based on high energy physics models, 
one finds that this uneventful story could be altered in several ways. It is important to identify 
these alternative possibilities and to relate them to specific observational signatures, that could be 
looked for in future cosmological data. 

An example of such alterations is the possibility that the coupling of the inflaton to other 
fields might affect its evolution. This is in particular the case if some field coupled to the inflaton 
becomes suddenly light at some special point along the inflationary trajectory. Indeed, whenever 
the effective mass of a field becomes small, there may be an associated burst of particle production, 
as thoroughly studied in preheating [l] (see also [2, [J) where this occurs while the inflaton is 
oscillating at the bottom of its potential. 

nn 

The same phenomenon can also take place during inflation, as pointed out in [4y-[6|], and such an 
event can temporarily trap or at least slow-down the inflaton by draining its kinetic energy. If such 
a trapping event occurs while observable modes are exiting the horizon during inflation this may 
lead to features in the power-spectrum (4 ] , as well as localised non-Gaussian spikes in the CMB, 
(see 0] for a recent comprehensive paper on this topic). By contrast, if this event occurs well 
after all the observable modes have exited the Hubble radius there might not be any observable 
signature (except perhaps primordial black hole formation on very small scales) if the trapping 
effect is exactly the same in all parts of our accessible universe. However, if the strength of the 
trapping fluctuates, this acts as an additional source of primordial adiabatic perturbations because 
the duration of inflation varies from place to place [jj. This modulation of the trapping can be the 
consequence of an explicit dependence of the coupling on another light scalar field as considered 
in |g]. 

In the present work, we primarily focus on the situation where the modulation is due to the 
fluctuations of the initial conditions of a multi-dimensional inflaton. Indeed, if one assumes that 
particle production occurs near a specific point in the multi-dimensional inflaton field space, then 
the trapping effect is stronger if the inflationary trajectory comes closer to this point during the 
grazing encounter. A slight change in the initial position of the inflaton thus affects the whole 
duration of inflation, depending on how much particle production slows down the inflatons and/or 
alters its initial trajectory. 
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The location in field space where additional degrees of freedom become light are often as- 



sociated with enhanced symmetries, which some authors consider beautiful 



These extra 



species/symmetry locations (ESL, ESP if it is a point) are a common occurrence in moduli spaces 
originating in string theory, 

seeJlO-Ql 

for a small selection of examples, and moduli fields can get 

trapped near such loci [9, 



1914231] (the string Higgs effect 



24J). If these fie 



and if several ESPs/ESLs are encountered, trapped inflation can result [91. I25I427I] (see also [281433 1 



ds drive inflation 



for related work); however, in this paper, we focus on a single grazing ESP encounter. Since infla- 
tionary trajectories are slowed down and derail from their prior course as the ESP is passed, one 
might say that beauty is distractive in our case. 

A modulated trapping event can easily produce significant non-Gaussianity of the local type, 
as stressed in . The current observational bound on the local type of the bispectrum from seven 
years of WMAP data is —10 < /nl < 74 at the 2a level [34|. The bounds on the trispectrum, which 



35(], and the CMB [3j438|]. 



are roughly tnl ^ 10 5 and \qnl\ ^ 10 6 , come from large scale structure 
If there is no detection, constraints with Planck should be tightened by an order of magnitude. 
Future CMB constraints on the trispectrum parameters were considered in 3S|. 



We find that a single close ESP encounter during inflation, as discussed in Sec. HH can yield 
the dominant contribution to the power-spectrum, see Sec. IIIIl if the coupling is of order one (we 
also discuss mixed scenarios, where S = V new jV total is neither close to zero nor one); this new 
contribution is accompanied by observably large non-Gaussianities. If the encounter occurs during 
the time-frame when large-scale modes are leaving the horizon, an additional bump-like feature 

1 Q, 0-4 



in the power-spectrum results ruling out the model for g 



42]; 

however, a subsequent 



encounter can successfully modulate large-scale fluctuations; the still present bump on smaller 
scales can, if observed in future experiments, be used to identify the time of the encounter and 
offer a consistency check. If non-Gaussianities are not observed, we can rule out the majority of 
ESP encounters, that is, we can severely constrain the interactions of the inflatons during inflation. 

In models considered in this paper, the dominant contribution to observables originates from 
the slowing down of the inflatons, Sec. Ill CI whereas the bending of the trajectory can be ne- 
glected, Sec. Ill Dl We derive approximate analytic expressions for the power-spectrum and non- 
Gaussianities, which we compare with a full numerical solution in Sec. Ill El We find that the 
analytic expressions are excellent for intermediate impact parameters fi ~ 10~ 3 Mp, which are 
needed anyhow if the additional contribution to the power-spectrum satisfies the COBE bound. 
For smaller impact parameters temporary trapping events can take place if the coupling is strong 
enough; however, since trajectories become chaotic and thus extremely sensitive to initial condi- 
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tions, we can rule out such encounters. A summary of the different regimes that we discuss is 
provided in Table HI 

We relate the likelihood of a single ESP encounter to the average inter ESP separation in (and 
the dimensionality of) field space in Sec. IIIIB1 where we also comment on the relationship to 
trapped inflation. Modest ESP densities with an average inter-ESP distance of y ~ 0.2Mp are 
sufficient in a two dimensional field space to render a grazing ESP encounter a likely occurrence 
during the last sixty efoldings of inflation. Thus, no particular tuning is required. 

A direct modulation of the coupling constant is considered in Sec. IIVI We compute the same 
observables in several cases, Sec. IIV A\ and comment on the possibility to use the modulaton as a 
curvaton in Sec. IIV Bl The two proposals, a grazing ESP encounter and modulated trapping, are 
closely related, but neither is a subset of the other. This is evident if the modulaton is promoted 
to a second inflaton, SecJV] In all cases we find comparable non-Gaussianities, and we provide 
consistency relations between the various observables n s ,r, fNL,nf NL ,TNL, 9nl throughout this 
article. 

Some technical aspects related to the validity of the approximations used to describe particle 
production (non-adiabaticity, sudden event) at an ESP are provided in the appendix. 

II. A GRAZING ESP ENCOUNTER 

We are interested in the effect that an ESP (Extra Species/Symmetry Point [j], [ijj]) encounter 
has on the trajectory, and subsequently on observables, in a 2D field space during inflation. An 
application of our results to higher dimensions is straightforward as long as only a single ESP 
encounter takes place. Throughout this section we set Mp = 1/ y (8irG) = 1. 

A. Particle production in a simple multi-field model 

Consider a simple multi- field inflationary model with two inflatons (fi, i = 1,2, canonical kinetic 
terms and a quadratic potential 




(1) 



In addition, we assume a coupling to another scalar field \ via the Langrangian 



int 



1 2 2/ ^ \2 
-j9 X VP - Vesp) • 



(2) 
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FIG. 1: Schematic of a grazing Extra Species Point (ESP) encounter with impact parameter /j in a two 
dimensional held space. If the dimensionality is bigger than two, the depicted plane (and coordinates Lp±, Lp\\) 
is defined by the location of the ESP and the tangent to the unperturbed trajectory at (p\\ = tp p . Note that 
ip± is not an instantaneous isocurvature field, but defined via a global field redefinition. In case of a simple 
potential, as in ([1]), trajectories are straight lines in the absence of an ESP. N — 60 cfoldings of inflation 
take place between \ip\ = </?* and \<p\ = </?end- 



For g > 0.001 radiative corrections to the effective potential are possible; these can be absent in 



super-symmetric theories where bosonic and fermionic contributions tend to cancel each other out 
[l|, Q], or in models within string theory 9, H . In the following we assume the absence of 
such corrections. When the inflationary trajectory comes close to the ESP, the \ particles become 
light 



m e J f =g\0-0 ESP \ + m* are 



(3) 



(in the following we ignore the presence of a small bare mass for the \ particles). Thus, if the 
impact parameter 



H = mm(\0(t) - <^ ES p|) 

is small enough (see Figures 1] and El) , particle production can occur, similar to preheating 



see [43|, 



44i | for reviews and 



(4) 



4514471] for multi-field preheating. For small impact parameters this 
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event is fast compared to the Hubble time during inflation, At p <C H, and can be treated as an 
instantaneous event in the following. 



The particle occupation number increase can be approximated by [l, 9| (see Appendix 1X1 for 
technical details) 

g 3/2 vf, /2 -nmt 

n ^=^r e vp ' (5) 

where v p = \<p(t p )\. The subscript p denotes the time of particle production when the distance to 



the ESP is minimal. Equation ([5]) is an analytic approximation valid for g > H p /v p [9] (always 
satisfied in this paper). The number of particles is diluted by the subsequent expansion 

n^n p {^y&(t-t p ), (6) 

so that backreaction is confined to a few efoldings at most (0 is the Heaviside distribution). 

Particle production is stronger the smaller the impact parameter /i is and the faster the trajec- 
tory is traversed. The velocity at which particle production occurs is easily found in the slow-roll 
regime when the slow-roll parameters e, = (V^/V) 2 /2 and rji = Vu/V are small (we use the 
shorthand-notation V t i = dV/ dipi ) . Given the potential in (pQ) the total slow-roll parameters along 
the slow-roll trajectory are simply eu = f)\\ = 2/(p 2 and we can solve the approximate Friedmann 
and Klein-Gordon equations to get the speed 

[2 

Vp ~ V 3 ?n = cons * ' ^ 
In our scenario, what matters is not the particle production in itself, but the influence of the 
particle production onto the motion of the inflaton. We thus need to compute the backreaction 
of these particles onto the inflatons. This is a complicated problem in general, but, in our case, 
it is reasonable to use the Hartree approximation, so that the equation of motion of the inflaton 
becomes (see Appendix ([B]) for details) 

dV dp Y , . 

difi dfi 

where 



p x ^n p [^y m ^, (9) 

is the energy density of the x field. This energy density originates from infusing kinetic energy of 
the inflatons into p x , which is usually small 1 , p x {t p ) <C = v p /2. 



1 The decrease in kinetic energy is accounted for in all numerical studies in this paper. 
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FIG. 2: Schematic of the change in trajectory of a grazing ESP encounter (no trapping, see Fig. [3] (c) for 
an example of a temporarily trapped trajectory): in the absence of an ESP the trajectory is aligned with 
the (f\\ axis and the dotted line is the bent trajectory due to backreaction of x-particles onto the inflaton 
fields. We approximate the production of %-particles as a single, instantaneous event at the point of closest 
approach. The amount of bending (|54[) and slowing down (|32[) depends sensitively on the impact parameter 
fi, which impacts correlation functions. Actual trajectories are plotted in Fig. [3] (a). 

B. Nonlinear curvature perturbation 

As we will quantify in the following subsections, any change in the trajectory caused by the 
ESP encounter, be it a change in velocity or an actual deviation of the trajectory, leads to a 
change in the number of efoldings, AA^. This change is sensitive to the initial conditions around 
sixty efoldings before the end of inflation at £*, which in turn causes additional contributions 
to correlation functions of the curvature perturbation on uniform density surfaces £, most easily 
recovered in the (non- linear) 5N- formalism [2,0]; these fluctuations are related to the fluctuations 
in the number of efoldings £ = SN, where 



where N SR corresponds to the number of efoldings along the slow-roll trajectory in absence of an 
ESP, while AN represents the change in the number of efoldings due to the ESP encounter. 

In this section we provide the relevant expressions based on the (W-formalism for correlation 
functions, before computing AA" numerically and, in certain regimes, analytically. 




(10) 



In our case, the total number of efoldings can be written as 



N = N SR + AN, 



(11) 



1. The power- spectrum 



Since the masses of <pi are smaller than the Hubble parameter, m <C H, all ifi carry nearly scale 
invariant fluctuations on super Hubble scales with power-spectrum 

.2 / rr \ 2 



(V kl V k2 ) = (2n)% j S^(k 1 + k 2 )P(h), P(k) = ^-P(fc), P = (^J 



(12) 



(In this paper ip± = <p l ). The index * denotes Hubble crossing of the pivot scale (&60 = a*-ff*) at 
N = 60 efoldings before the end of inflation. 

Expanding the number of efoldings at first order as a function of the inflaton perturbations, 
6(fi, or alternatively 8<p\\ and 5ip± (where (pu and (p ± are the components respectively parallel and 
orthogonal to the trajectory without particle production; see Fig. [T|), the curvature perturbation 
can be written as 



( = 5N = 5N SR + 5AN = N\\6<p\\ + AN :± 5ip± , (13) 



with the notation 



dN dAN 

= AN ^i^- (I4) 

(we employ corresponding expressions for higher-order derivatives). Here we used that due to the 
geometry of the ESP encounter, the impact parameter \x and consequently A./V are only sensitive 
to changes of the initial field values perpendicular to the trajectory. Note that the perpendicular 
component of the inflaton plays a similar role as the modulaton field introduced in [8(. 

Hence we can write the power-spectrum as a sum of two contributions, the usual slow-roll 
contribution and a new one caused by the ESP encounter, 



— + AN 2 

2e* 



(15) 



,2vr J 

where e* = —H/H 2 is the usual slow-roll parameter. Here AiV is comprised of two contributions 
that decouple at first order: one due to the slowing-down (subscript "si") of the fields and the 
other due to the geometric change of the trajectory (subscript "ge"), 

AN = AN sl + AN ge . (16) 

Defining 

S = (17) 



the scalar spectral index is given by [5l|] 



1 55 = " (6 " 4HK + 2(1 " + 2Evl - (18) 



where r/jj and r]\ correspond to the usual slow-roll parameters parallel and perpendicular to the 
trajectory evaluated sixty efoldings before the end of inflation. For 3 — > 0, the slow roll resu 
recovered, while for larger 3 the model is similar to the mixed inflaton/curvaton scenario [52 



t is 
or 



modulated traping [8|] 2 , see forward to Sec. IIII1 The tensor-to-scalar ratio can be computed to 



5l| 



r = lfe*(l-E). (19) 

2. The hi- and trispectrum 

One can use the 5iV-formalism to compute the non-Gaussianities of the curvature perturbation 
associated with the trapping effect. Let us first recall the definition of the local non-linearity 
parameters Jnl (bispectrum) as well as t^l arid g^L (trispectrum). For a multi- field system, the 
expansion of the number of efoldings in terms of the scalar field fluctuations can be written as 

( = SN = N,i dip' + ^N }ij V V + \N, ijk dip' V Sip k + . . . , (20) 

where we use the implicit summation convention for the field indices k, . . . and the notation 
Ni = dN/dip%, N,ij = 9 2 W/VVk etc. 

If all of the scalar field fluctuations are Gaussian, with the power-spectrum (|12p . one finds, using 
Wick's theorem, that the bispectrum, i.e. the three-point function in Fourier space, is given by 



(c kl Ck 2 c k;i ) = (2^) 3 <5( 3 ) n£kiW(ki,k 2 ,k3), 



6 

£? c (ki,k 2 ,k 3 ) = ^hh[P(k 1 )P(k 2 ) + P(k 2 )P(k 3 ) + P(k 3 )P(k 1 )], (21) 



with the non-linearity parameter 



6 _ N ti N tj N# 

5 /NL " (N, k N*)* ■ (22) 



Similarly, the trispectrum, defined as 



(C kl Ck 2 Ck 3 Ck 4 ) c = (2^) 3 5 (3) (l>^ 7c( k i' k 2,k 3 ,k 4 ) , 



(23) 



In this reference, r/1 <C was assumed. For us rf^ = due to JTJ. 
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can be written in the form 53] 



54 

r c (ki,k 2 ,k 3 ,k 4 ) = r NL [P(k 13 )P(k 3 )P(k 4 ) + 11 perms] + — g NL [P{k 2 )P{k 3 )P(k 4 ) + 3 perms] , 

(24) 

with 

NijN^N'iNk 25 N^N^N* N> k 

TNL " (NjNJf ' 9NL ~ 54 (NjNJf (25) 

and where k\ 3 = ||ki + ka||. 

In the present case, one finds that the non-linearity parameters due to the modulated trapping 
are given by 

5 AN* ± AN ±± 5^ 2 AN ±± 
JNL 6(AiV 2 +M) 2 6" (AiV ± ) 2 ' 1 J 

' ill 

AN%AN% ± (AN ±± r ^ _ 36 ^ 2 
^ (AiV^ + iV 2 ) 3 (AiV i± ) 4 " 25^ Jnl ' 1 ' 

25 AJVf_ L AJVj.j. i . 25 AiV ±±± „ 3 
57VL 54 (AJV* + JVp 3 54 (AiV ± ) 3 ^' 1 j 

The inequality that t^l > (G/nl/^) 2 , which here follows from the fact that 3 < 1, is true much 
more generally J54]. However it might be violated in special models in which the "loop" corrections 
are very large [551 ] . 



C. Slowing down 

Our next goal is to estimate AN s i analytically, closely following 8| . Assume that after the ESP 
encounter, the change in ip j_ remains small compared to cpu and (i and consider that |(/?|| — ip p \ S> /U. 
In this case, we can approximate the equation of motion (jHJ) for ipu by (see (1B5|) in the appendix) 

<p { \ + 3H p <p {l + = gn p Q(t - t p ) . (29) 

We would like to solve this equation over a brief time interval of a few Hubble times only; in this 
case we can treat H « H p = const and integrate (f29|) to 

¥>|| « -v p + gnpe" 311 ^-^^ - t p ) , (30) 

were we imposed \cp^ \ = v p at t p . By denoting 

A n (t) = n (t,g^0)- n (t,g = 0), (31) 
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we get 

Ap|| =gn p e- 3H ^ t - t ^(t-t p ), (32) 

which gives the amount by which the field is slowed down after the ESP encounter. Thus, the 
trapping event delays the field by 

A n = f A n dt = %wi> (33) 

and the net change in the number of efoldings becomes 

AN sl « (34) 

fl 5/2 1/2 M 2 

- 5 Vp r e~^. (35) 



9^ p (2vr)3 

Note the presence of the exponential, which is absent in the one dimensional case. 

A small change in the position of the inflaton at time t* induces a variation of the angle f3 
between the inflationary trajectory and the i^i-axis (see Fig. [1]), which entails a change of the 
impact parameter 

5fM ~ <p p 5/3 —(p* ± . (36) 
In this regime, the power-spectrum (fT5j) picks up a piece that is boosted by 

AiV^ « (-Y^-jfv^ 2 ^ - A^xe-i* , (37) 

where we have defined 

A _ g 5/2 ^ ^ c _ 2 vp m_ ^ x _ ->ml_ _ (38) 

if p 9(27r) 3 ' 99* Y v p v p 

Moreover, the additional contributions to the non-linearity parameters due to the slowing down 
effect become 

, 5 A 3 C 4 x (x - i) e~ 3x , s 

f S NL~ r , * 2j — 2 , (39) 

NL ~ ^T3 ' ( 40 ) 



A 2 C 2 xe~ 2x + f 



25 A A C & x 2 (x - |) e 



- l,r 



k (41) 
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Here, we also assumed that effects due to the geometric change of the trajectory are much smaller 
than the ones caused by slowing down, see Sec III Di 

Under which conditions can we trust these analytic expressions and when are they observable? 
The additional contribution to the power-spectrum is maximal at 



with a corresponding amplitude of 



^(^)-^f^4-1.2xlO-V, (43) 

where we used the slow-roll approximation ip p H*/((p*H p ) ~ 1 as well as (3 <C 1 (see Fig. [T]). 
Given a large enough coupling g, it is possible that this contribution dominates, saturating the 
COBE bound. To find the coupling beyond which the trapping effect can dominate we solve 
r s i/V SR = 1 to 

1/2 



9d 



(^ 2 ) 1 / 3 ((2vr) 5 e ) 1/6 (0 « 3.9(^)1/3. (44) 



If g > gd, there exists a range around \i c for which the amplitude is larger than the usual slow-roll 
contribution. Of course it is then this contribution which needs to be tuned, i.e. by carefully 
choosing ip* ± . We investigate this model, including a discussion of non-Gaussianities, in Sec. IIII1 

If the coupling becomes too large, a qualitative change in the trajectory occurs: instead of 
being merely slowed down, the velocity reverses and the trajectory gets temporarily trapped. This 
takes place when the force due to backreaction is bigger than the one stemming from the classical 
potential, 



dV 



dipi 



< 

p 



Px 
dipi 



„5/2 3/2 2 

g v p e —%- (45) 



where the subscript p refers to evaluating a quantity at the time of particle production t p . Solving 
for g in the \i = case yields 

/ /s\ 3 / 4 \ 2/5 

9trap = (m^)i/5 I (2^)3 y j w 3.2(m^)V5 , ( 46 ) 

which is bigger than the coupling at which the power-spectrum can be dominated by V s i 3 . 



3 9d ~ — — — ~ <Pp*' < 1 f° r H 3 * — 16 and m = 10 B ; since inflation ends for 1.4 ~ ip < ip p < tp* and the 



Strap 



inequality becomes stronger for smaller m we get <?tra P > gd- 
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However, even for large g > g^p, the slowing down effect is well approximated by the above 
analytic expression if the impact parameter is bigger than 



A* > A*an = 




(47) 



This parameter is usually bigger than the one for which the power-spectrum has its peak, but low 
enough to allow for the use of the analytic approximation in Sec. IIII1 



D. Geometric change of the trajectory 

In this section we would like to estimate the path change compared to the trajectory in the 
absence of an ESP. The trajectory bends towards the ESP after particle production took place. 
Our goal is to estimate the maximal deviation in the perpendicular direction A</?j_ under the 
assumption that A(^_|_ <C /i, which in turn leads to a slight enhancement of the number of efoldings 
AiVi, in addition to the slowing down effect discussed in the previous section. Since /i > Ac/?_|_ 
(see FigE]), the effect is small compared to the slowing down effect. It becomes much stronger for 
large coupling g > g trap and small impact parameters < ji < fi an , for which it is challenging to 
find analytic results. We compute such cases numerically in Sec. Ill El enabling a comparison of 
our analytic approximations to the full solution. 

There are two different regimes for <p± <C fi: in the first one, back-reaction is strong and the 
trajectory bends, which requires \ip\\ —(f p \ ~ v p (t—t p ) <C H, and we can approximate the equation of 
motion (JSj) by (for ease of notation, we put the ESP above the trajectory in this section (/2_i_ E sp = +A 4 
with fi > 0, opposite to Fig. [1] and Fig. [2] but in line with App.[Bl i.e. eqn. (|B6p ) 

(px_ + 3H p <p± + m 2 tp ± ps gripe' 311 ^-^ , (48) 

where we treat H p « const, since we are only interested in following ip± until v p (t — t p ) = [i (a 
fraction of an efolding). Ignoring the classical potential (m = 0), the above can be integrated to 

Hence, the excursion at t — t p = fi/v p is of order 

^-lmy 2 ^ = A v 2ip ^ (50) 

with n p from ([5]), we expanded the exponential for small v p (i/H ~ 2fi/ip p <C 1 and we used Ai^y 
from (|33p . We see that A(p± <C H for small fi. The corresponding velocity in the _L-direction is 
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approximately 



^± = y±(t)~ 2g ■ (51) 

Once — tp) > (J,, the source term on the right hand side of (jHJ) becomes unimportant quickly 
and the second phase commences: the field rolls freely until the classical potential and Hubble 
friction bring the trajectory back onto a slow-roll attractor. This takes place roughly within 
one Hubble time, see Fig. O panel a. The total excursion in the perpendicular direction can be 
estimated by solving 

tp± + 3H p p ± » (52) 

with (p™ 1 ~ (we checked that using tp± from f|50j) does not affect results significantly) and 
(p 1 ™ 1 = v± from above to 

^W-^(l-e- 3 ^-^) , (53) 

and taking the large-t limit (which maximises <f±(t) in the limit that m± = 0), we find 

3 

Aip± « A^n^p- . (54) 
Thus, the trajectory is prolonged by roughly 4 

3 

As ~ Aip± « A<^||/x^p- , (55) 

and the corresponding change in the number of efoldings can be estimated, up to factors of order 
unity, to 

AN „ EAl . (56) 
v p 

For small impact parameters, this change in the number of efoldings is much smaller than the 
corresponding one from slowing down, 

AN 9 e _ . 3^ P M _ 



AN, 



si 



4 An alternative viewpoint might be instructive: during the first phase, the trajectory is bent by the angle tan(a) = 
¥>±{t)/M>, with respect to the unperturbed trajectory; since this angle is small we have a ~ A(p\\fi(fp9/8. Since 
the overall speed decreases in a grazing ESP encounter and the decrease is not strong in the cases of interest, see 
Fig. [3] panel b, we take the slow-roll speed v p as an approximation for the speed along the intermediate trajectory 
before a slow-roll attractor is reached again (this provides an upper bound on effects). Then the difference in 
path-length during the excursion is of order As ~ v p /H p — I, where I — v p cos(a)/ H p ss (1 — a)v p /H p . Plugging 
in a gives ~ jfi<p P , consistent (up to factors of order one) with the above estimate. 
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Hence, as long as \i <C 2/(</? p 3), slowing down effects are dominant. Furthermore, since the ratio 
of the derivatives is 

A^^--^^- 1 ) ' (58) 
the slope is dominated by the slowing down effects for fj, an < fj, <C 2/(3<p p ) where fi an = few x 



\Zv p /(irg) (the same holds true for higher-order derivatives up to factors of order unity). For 
small fj, and small g the geometric effects can be bigger than slowing down effects at the level of 
correlation functions; however, in this regime the usual slow-roll contributions tend to dominate, 
since corrections scale as fi. This is the reason why, for \i !> f^ani 

we are justified to ignore AN ge 

in Sec. Ill CI when computing the power-spectrum and the non-linearity parameters. Then, the 
additional contribution to the power-spectrum and higher-order correlation functions caused by 
the geometric change in the trajectory are of order 

-P 9 e ~ AWe-^fx-jV, (59) 

(60) 



2vr, 

' NL '~ fi / 

b (A 2 C 2 xe~ 2x + f 

t nl ; ' ( 61 ) 

A 2 C 2 xe -2x + K 



25 A^B^e-^ (x 2 - 3a+| ) (x - ±) 3 
^A 2 C 2 xe~ 2 - + 



* Z ~ . (62) 



where we defined 



B=^M. (63) 

If g < g trap , so that the analytic approximations are valid for all and fi <C /i an , so that 
|A./V*[| <C \AN 9 ?\, one only needs to replace the denominators above by appropriate powers of 



(C 2 A 2 B 2 e~ 2x (x - 1/2) 2 + (pi/ 4). Note that AN g [ changes sign for (i < y / v/(2-Kg), indicating that 
increasing the impact parameter from zero up first leads to an increase in the number of efoldings, 
in accord with the trajectory being bent, before the weakening of the slowing down effect decreases 



Ntotai and, for fj, > y/v/(2irg), trajectories show less of an excursion again, see Fig. [3] (a). 



E. Discussion and numerical studies 



We encountered several different regimes with qualitatively different behaviour, see table [J for 
a summary. All regimes can be treated numerically without invoking the slow-roll approximation 

16 



Coupling Constant Impact Parameter 


Description 


9 < Scobs in (|64 


1 A 4 = any 


rCOBE > PbSP 


g d < g in (gU) 


/x near fi c in (glj) 


Pinfl < V ESP 


.9 < ffmax in dZ3) 


/x = p in ([75)) 


The non-adiabaticity parameter in (|A2|) is bigger than one 
and particle production is fast (this paper). 


3 > Smax 


/i = /l in (f?3"|) 


The approximation in App. C of [9] could be used, 

but NG are expected to be too large for m ~ TO CO BE/few. 


gd< g < Strap in 


0BJ n < 2/(^3) 


Slowing down effects dominate over geometric effects. 




/i = /i 


If m < m CO BE in (l66l) then 7-cobe is saturated by V s i at p. 




/i > 2/(^3) 


Geometric effects dominate over slowing down effects, 

but they are usually negligible compared to slow-roll effects. 


Strap < s 


/x < fx an in (jUJ 


Temporary trapping can occur, 

N(n) can become non-analytic, Fig. |4](d). 




Mon < fJ- <^ 2/(ip p 3) No trapping, slowing down effects dominate. 




n = n 


If m < to CO be then 7-cobe is saturated by V s i at p,. 




H > 2/(^3) 


Geometric effects dominate over slowing down effects, 

but they are usually negligible compared to slow-roll effects. 


0.8 < g < 2.5 


/i = /i 


'Pcobe is matched by V s i in Sec. Mil and non-Gaussianitics 
are observable (we use 0.02 < m/m C oBE < 0.2). 



TABLE I: Different regimes depending on the coupling constant g and the impact parameter ji. 



(which can be violated during the ESP encounter) or the small p approximation, as long as we 
assume a single, instantaneous particle production event. The latter can be a good approximation 
even for temporarily trapped trajectories, since the velocity decreases after the first encounter with 
the ESP, v < v p , and particle production is less efficient subsequently, ([5]) (see Fig. [3] (c) and (d) 
for a concrete case). 

First, we would like to investigate some trajectories in field space and how they compare to our 
expectation. In Fig. [3] (a) we can see the that the bending of a trajectory grows as p is increased 
from zero, before they remain essentially straight for larger impact parameters. This is in line with 
our predictions in Sec. Ill PL where we observed a sign change in the second derivative of AN ge in 
(1601) . The corresponding velocity change is plotted in panel (b). At the particle production event, 
only a small fraction of the energy is transferred to x particles, less than 1% in Fig. [3] (b), but 
the velocity decreases subsequently due to backreaction of x particles. As the universe expands, 
p x , and thus the force due to backreaction, redshifts as a~ 3 , causing the field to speed up again 
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FIG. 3: (a) and (b) show the trajectory and corresponding velocity v = V01+02 ( over t/tjj ~ trrup p / 'V6) 
directly after an ESP encounter (x particles are produced instantaneously at t = only) with g = 1 < 
5tra P ~ 1.8 for a range of impact parameters \i (the same as in Fig. [3] (b)). The field slows down while the 
trajectory gets slightly disturbed, panel (a). The slowing down effect, which dominates correlation functions, 
is monotonically decreasing with increasing fi, and v approaches the slow-roll speed again after a few Hubble 
times, panel (b). (c) and (d) show a temporarily trapped trajectory with corresponding velocity, with an 
increased coupling g = 3.2 > g tlap . The fields oscillate temporarily around a non-parabolic minimum of the 
effective potential, and the resulting trajectories are highly sensitive to initial conditions. This leads to a 
non-analytic dependence of the number of efoldings on initial conditions for small impact parameters, Fig. |4] 
(d). In all plots m = 0.2 x tocobe, defined in (larger m increase SWp)- 



and approach the slow-roll speed v p in ([7]) after a few Hubble times. The maximal decrease of 
the velocity is monotonically decreasing as \x increases, since the exponential suppression in (135|) 
becomes important and as a consequence, ffj L in (|39p is always positive. 

The coupling in panels (a) and (b) of Fig. [3] satisfies g < g tra , P , preventing a trapping of tra- 
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FIG. 4: The total number of efoldings N in the presence of an ESP divided by the number of efoldings 
in the absence of on ESP, N noESP , is plotted over the impact parameter [i, for (a) g = 0.5, m = to CO be 
(slow-roll contributions dominate the power-spectrum), (b) g = 1, m = 0.2 x m CO BE (the slowing down 
effect can dominate the power-spectrum), (c) and (d) g = 3.2, m = 0.2 x tocobe (for small [i the trajectory 
becomes temporarily trapped, the system becomes chaotic and N(fi) non-analytic, see panel (d); for (i > /i n, 
trapping does not occur and the slowing down effect dominates again). 

jectories. We chose g = 1 as a representative coupling, since it is the most natural value for a 
dimensionless parameter; it is large enough, g > ^cobe (see forward ()64|) ) and ga, for the power- 
spectrum to be dominated by the slowing down effect, Fig. H] (b) and Figj5] (a). If the coupling is 
increased above g trap ~ 1.8, trajectories can get temporarily trapped for small impact parameters 
fJ- < ^an- An example of such a trajectory is shown in Fig. [3] (c) with the corresponding velocity 
in (d). 
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FIG. 5: The additional contribution to the power-spectrum normalized by "Pcobe as well as the non- 
linearity parameters /nl, tnl and Qnl are computed numerically (crosses) and compared to the analytic 
approximations (red line) in Sec. Ill Cl which lead to (|57|) and ([59|) - (|1T|) using 3 = 1 and m = 0.2 x m CO BE; 
slow-roll contributions are negligible and as a consequence, the impact parameter \x needs to be tuned in 
order to satisfy the COBE bound (the normalized power-spectrum in (a) needs to equal one, singling out 
two values of fj,), see forward to Sec. IIIII The discrepancy between numeric and analytic results for large 
and small fi are due to additional effects caused by the geometric change of the trajectory, see Sec. Ill Dl 
which arc negligible in Sec. IIIII where V jVcow = 1 at the larger of the two possible impact parameters. 
Rounding errors amplify quickly once "P/Pcobe "C 1 (negligible in Sec. IIII[) . 

If trajectories get trapped 5 , they oscillate around the minimum of the effective potential close 
to the ESP. These trajectories are highly sensitive to initial conditions, as expected for a two 

5 The way particle production is described as an instantaneous event here and elsewhere [sj is not a good approx- 
imation to describe trajectories during the actual trapping. Thus, plots like Fig. [3] (c) and (d) as well as Fig. |4] 
(c) and (d) should be taken with caution. However, the statement that trajectories get trapped and that they are 
sensitive to initial conditions is robust. 
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dimensional non-harmonic oscillator (the system is chaotic; see [23(] for related work). For example, 
a trajectory can change qualitatively from having three revolutions to two revolutions at a critical 
impact parameter. This leads to discrete jumps in the number of efoldings, rendering N(fi) non- 
differentiable. For such a chaotic system the 5N formalism can not be applied. This is evident in 
Fig. 2] (c) and (d): for large fi, trajectories are not trapped and the shape of N(fx) is similar to the 
ones in panel (a) and (b) where g < g tT a. P , but for small [i the curve in panel (c) becomes irregular. 
A zoomed in version of this small ji range in panel (d) reveals that N(/i) is not smooth at all. An 
investigation of chaos in such a system is interesting 6 , but goes beyond the scope of this article. In 
the following we focus on cases without temporary trapping > fj, an if g > g trap ). A remnant of 
this sensitivity can be observed in regions where geometric effects become important, see forward 
to Fig. [5] (d), where a discontinuity in the second derivative of N causes a divergence in g^i for 
small //. 

If slowing down effects dominate over geometric effects 7 , \x an < [l <C 2/(ip p 3), N(/j,) is smooth. 
One can then tune the impact parameter to match the COBE normalization and compute corre- 
sponding non-Gaussianities, see Sec. IHII and Fig. [6] In order to make analytic predictions, we focus 
entirely on slowing down effects as described in Sec. IIIII This is an excellent approximation, which 
can be seen in Fig. [5] (a)- (d), where the full numerical solution (crosses) for the power-spectrum and 
the non-linearity parameters /nl, tnl and gNL are compared with the analytic approximations 
(f3T|) and ([39|) - (|1T|) for g = 1 and m = 0.2 x m COBE . We chose m < m COB E and the most natural 
coupling g = 1, since we focus on such cases in Sec. IIIII The differences for small and large \i are 
caused by geometric effects becoming important. However, it is clear that around the maximum 
of V/V cobe the analytic approximations are excellent. 



6 For example, do intermediate regions of smooth N(n) exist? What is the fractal dimension of the non-smooth 
set? Etc. Note that in order to answer these questions, one needs to track subsequent particle production and 
the amplification of numerical errors. All numerics in this paper were done with MAPLE9's standard integration 
routines, 14 digits accuracy (we varied this accuracy to check that no conclusions in this article are sensitive to it) 
and a single particle production event. 

7 In the case of trapped trajectories, such a separation does not make sense anyhow, since neither of the two analytic 
descriptions in III CI or III Dl are applicable. 
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III. CURVATURE FLUCTUATIONS FROM A GRAZING ESP ENCOUNTER 



A. Confronting the model with observations 

If the coupling of the extra species particles is strong, g > g<i in (144j) . the power-spectrum can 
be dominated by fluctuations in tp± due to the slowing down effect (Sec. Ill CP provided that the 
impact parameter is close to fi c in (|42|) . but less than \i <C 2/(ip p 3) at which point the geometric 
effect of Sec. Ill Dl becomes important, according to (|58|) . To satisfy the COBE normalization via 
the slowing down effect alone, we need 

9 > Scobe = (e(2^) 7 3 4 P COB E) 1/6 « 0.77 > g d , (64) 

where <? C obe is defined such that V s i(fJ, c ) = "Pcobe- For g > gj ~ 0(1) and tp p ~ 10, the inequality 
fj, <C 2/(ip p 3) implies fi <C 0.1, which is easy to satisfy. Since g can be larger than <? trap in (j4"6j) . and 
we would like to focus on trajectories that don't become temporarily trapped (then AN± can be 
non-analytic, Sec. Ill E|) . we require fx > [i an from (|4T|) in this section. Under these assumptions, 
the slowing down effect is recovered well by the approximations in Sec. Ill CI 

Since we do not want the COBE normalization to be saturated by slow-roll contributions, that 
is by perturbations in cp» , we require |34| 



T- « ^cobe ~ 2.41 x HT 9 , (65) 
or 

m < wicobe = 6.2 x 10" 6 . (66) 

This amounts to less fine tuning of the potential, since only an upper limit is imposed. On the 
other hand, the impact parameter needs to be just right to saturate the COBE bound, that is 
x = irgfi 2 /v p needs to solve 



xe ~ 2x = D , (67) 



where we used (1371) and defined 



D - 1|) 2 -^)^ (68) 

^(f) 2 ^ (09) 

2 6 3 4 7T 7 - 

S ^PcOBE ~ 1-6 X 10 7 ^P CO BE (70) 
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with A and C from (|38p . If slow-roll contributions are negligible S — > 1 (the proper value for E 
is kept in all plots). For g > <7cobe, D is smaller than the maximum value of xe~ 2x at x c = 1/2 
(corresponding to /i = /i c ), that is D < l/(2e); then equation (|67|) has the two real solutions 

Xi = -hVi{-2D), (71) 

for i = 0, —1, where Wj is the Lambert function, i = leads to an impact parameter smaller than 
fi c , and % = — 1 to a bigger one; we focus on i = —1 primarily because impact parameters smaller 
than n c lead to a large negative /nl (see Fig. [5]), that is observationally disfavoured 8 and define 

x = x_i = -iw_i(-2I>) . (72) 

This lower branch of the Lambert function is defined for arguments in the interval [— 1/e, 0], but 
it has no simple series expansion at the boundaries (we are particularly interested in arguments 
close to the lower boundary of the interval). 

How well do we have to tune x and thus the impact parameter to be within 1% of 
P COBE ? From Eq. ([6?]) we find xe~ 2x = D x (1 ± 1CT 2 ) leading to x = x ± 5x with 
Sx = \10~ 2 W^ 1 (-2D)/(2D(1 + W-i(-2D)))\ « 0.08 ~ CflO" 1 ), where we performed a Taylor ex- 
pansion around x and took 9 L> ~ 10 _1 so that x ~ 1.3 in the last step; this corresponds to a tuning 
of about 6% of x, or about 3% of 

|U = |U c \/2^. (73) 

We observe /z > /i c 3> 5ipJ_ M as long as (p py fg <C > 2/- v /m COB E ~ 800. For large m close 

to m CO BE and ~ 10, the required tuning of the impact parameter 5ft becomes uncomfortably 
close to dtp®^ 1 if g ~ 1; hence, we require the inflaton mass to be at least a factor of 5 smaller than 
"Icobe to be on the safe side. 

If m < m COBE , then "Pcobe is matched by a mix of slow-roll and slowing-down contributions. 
In the limit m/m COB E *C 1 the slow-roll contributions become negligible, S — )■ 1, and the scalar 
spectral index (fTKj) becomes 

n s - 1 « (74) 
to first order in slow-roll parameters, since 77* = e*. An n s this close to one is observationally 



disfavoured 



34l | at the two sigma level . However, a scale dependence in line with observations 



8 We also want to prevent a tuning of /j, smaller than quantum fluctuations in any given Hubble time, Sip^ 1 ~ 
H/(2n). Requiring [i c /d(p c l M 2> 1 amounts to tfp^/g where m < mcoBE- 

9 In the limit D — > the tuning becomes more severe since Sx — > as well. 

10 For more discussion of the spectral index in a related scenario see [5(| • At second order in slow roll, n s — 1 = — 10e 2 /3 

m. 
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can easily be introduced by considering more general polynomials in (pQ) with e* ^ r]*j_ (we leave this 
to future studies). The tensor-to-scalar ratio in (|19|) approaches zero in the limit 3 — > 1 regardless 
of the potential. Thus, a detection of gravitational waves of order r ~ 0.1 would rule out a grazing 
ESP encounter as the dominant contribution of the power-spectrum. 

As discussed in appendix [A] the impact parameter needs to be smaller than /i max = yVp/g, to 
guarantee a large non-adiabaticity parameter and allow for particle production to be treated as 
instantaneous. Since /2//x max = l/v2~7r, we see that x < it has to hold. Solving x(D) = it using 
(|72p . leads to D min ~ 0.0059, which in turn leads to an upper bound on the coupling g via (|70p of 

/ 2 6 3 4 vr 7 \ 1/6 

9 < Smax = H— PcOBE » 1.4, (75) 

\ -^min / 

where we used 3 ~ 1 in the last step. Therefore results for coupling constants above g max and 
fx = fi, i.e. in Fig. [6j should be taken with caution, since our approximations start to become 
unreliable . 

To briefly summarize, given a value of the coupling in the interval 1.4 « g max > g > <? C obe ~ 0.77 
and m ~ m COBE /few, the impact parameter needs to be /i = p, within a few percent in order for the 
power-spectrum to be dominated by fluctuations in (p± while satisfying the COBE normalization, 
with n s — 1 ~ (different potentials can provide a red tilt) and unobservable small tensor-to-scalar 
ratio r. Next to the suppression of primordial gravitational waves, how can we tell this mechanism 
apart from a simple slow roll model? 

Since the above mechanism is contingent on a relatively strong dependence of the number of 
efoldings onto the isocurvature direction ip±, one should expect large non-Gaussianities 12 . This is 
indeed readily checked, by plugging x of (|7ip into the expressions for the non-linearity parameters 
(|39p -(j4ip. A plot of these parameters over g > 1.1 x ^cobe is given in Fig. [6j where we chose 
ip p ~ 11, so that the ESP is encountered around 30 efoldings before the end of inflation; we varied 
the inflaton mass in the range 0.02 < m/m COBE < 0.2. 

First, we observe that /jvl is generically observable (few < Jnl < 100) for <7 C obe < 9 < 0ma X 
and m = m COBE /few. fjy^ increases monotonically with increasing g and decreases as the inflaton 
mass is decreased. As a consequence, for low inflaton masses, a relatively broad range of coupling 

11 One could employ the iterative analytic approximation scheme in Appendix C of 0| (an expansion for small non- 
adiabaticity parameters) if g ^> <7 ma x and /j, = p,, but since we expect non-Gaussianities to grow uncomfortably 
large, we refrain from applying it here. 

12 Since we use a quadratic potential and equal masses for the two fields, the scale dependence of Jnl, see forward 
to Sec. IIVI is negligible. If a different power-law were used, a scale dependence of the non-linearity parameters 
would follow. 
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FIG. 6: The non-linearity parameters Jnl, t~nl and g^L in ()39[) - (j41|) from the slowing down effect, evaluated 
at fi = p, in (|73j) for which the COBE normalization is satisfied; slow-roll contributions are negligible for the 
power-spectrum, since m <C tocobe = 6.2 x 10 -6 ; we vary 0.02 < m/m cow < 0.2 and 1.1 x gcow < 9 < 2.5. 
For masses not too small compared to dicobe, non-Gaussianities are observable and keep increasing with g. 
Results are indicative even for g > .g t ra P as long as /2 > fi an (satisfied above), but since <7 t ra P > 5max in (|75|) . 
the sudden particle production approximation ceases to be valid and effects may be overestimated. 



constants permits observable non-Gaussianities that are not yet ruled out, but results become 
untrustworthy above g > (7 max . Values consistent with current observations are located around 
g ~ 1, the most natural value for the coupling constant, so no extra tuning is needed. The 
non-linearity parameters of the trispectrum are naturally of order f^ L (Fig. [6]): 
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tnl/Inl i s independent of g and only weakly dependent on m via H* in E, while qnl/ Inl 
also depends on g via x. It is clear that measuring both the bi- and trispectrum would help to 
distinguish this model from other scenarios 



54|, 



581 ] (see [59( for relationships between higher-order 



parameters). See the introduction for constraints on these parameters. 

We put the ESP encounter around N p ~ 30 efoldings before the end of inflation, well after 
observable scales left the horizon, for a good reason: if the ESP encounter were to take place close 
to or within the observational window, back-scattering of x particles onto the inflaton condensate 
would lead to an additional bump-like contribution to the power-spectrum and non-Gaussianities 
via IR cascading {tII^I-I^. This unobserved contribution would dominate the power-spectrum (and 
non-Gaussianities) for g > 0.1 if /i is small 13 , well before the ESP encounter starts to matter. But 
as long as the ESP is grazed later, the bump like feature is on smaller scales which are not tightly 
constrained by observations (but might be in the future). Choosing N p different from 30 leads to 
changes of order 0(1) in the constants A, B, C in (|38p and (|63p . but x in (|72p is insensitive, since D 
in (j68|) depends only weakly on tp p . As a consequence, the non-linearity parameters in (|39p - (|4ip are 
simply rescaled in the S — > 1 limit (Jnl oc A oc 1 / 4> p oc 1/ y/~Np~ and tnl, 9nl oc A 2 oc 1/0? oc l/N p ) 



j p w. J- 1 -y Jip 

and the consistency relations in (|76p and (|77p remain unaffected. 

The observation of a single bump caused by backscattering in the vicinity of k^ = ^Jgv p exp(Nb) 
where Nj, = N — N p , approximately 14 given by 

/7re\ 3/2 k 3 -f 4 t 9 ^ 2 
n-A 6 (f) F e 2 ^e , (78) 

6 



with g ~ 1 and ~ lO -6 ^ 15 / 4 |40j, |4l|], is a smoking gun of a grazing ESP encounter. Note that 
ngji 2 /v p = x = —W-i(—2D)/2. The location of this feature determines when the ESP has been 
encountered, and the height offers a consistency check of the model, since no free parameters are 
left. This bump is accompanied by large, potentially detectable, non-Gaussianities Q, Q] in the 
vicinity of kf,. 

It should be noted that all of the above arguments remain valid in higher dimensional field 
spaces and more general polynomial potentials (as long as all fields have the same potential); then 
the two fields that we consider above parametrize the plane spanned by the unperturbed trajectory 
and the ESP, while contributions from additional perpendicular directions are ignored. 



13 The bound g < 0.1 is derived in the fi = limit. As fj, grows, the exponential suppression in ((78} suppresses the 
amplitude so that larger oc g 15 ^ 4 and thus g are possible. In our case, ngp^/vp = x = —W-i(—2D)/2 is close 
to one so that the bound is not raised much (for example, x ~ 1.3 for D ~ 10~ so that g > 0.14 instead of 0.1). 

14 This approximation does not retain oscillations in the tail of the exponential [iol. l4lt|. 
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If the potential is more complicated, one may still compute additional contributions to correla- 
tion functions in a similar manner, with the appropriate changes in i.e. the relationship between \x v 
and ip*j_ in (|36p. If the trajectory goes through an effective single-field stage after the observational 
window is passed but before the ESP is encountered, all sensitivity to the initial field values are 
erased and no additional super-horizon fluctuations are imprinted onto cosmological fluctuations 
when the ESP is passed. By single-field stage we mean a regime where all trajectories are funnelled 
through a valley in the potential such that all but one direction have a heavy mass. 

So far, we just assumed that an ESP is closely grazed, but how likely is such a single encounter? 



B. Several ESP encounters and connection to trapped inflation 

If ESPs are dense we expect more than one close encounter during the last sixty efoldings 
of inflation. If x-particles dilute sufficiently in between encounters and none of them involves 
actual trapping, the contributions to the power-spectrum and higher-order correlation functions are 
additive, and the results of the previous section can be used directly. However, if encounters start 
to overlap, that is if one or more ESPs are passed per Hubble time, their combined backreaction can 



lead to a reduction of the speed throughout inflation, enab 
phenomenon is sometimes referred to as trapped inflation 
in one dimensional 
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ing steeper potentials to be used. This 



mgsteec 
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271 ] . which has been investigated 



261 ] field spaces. 



271 ] and higher dimensional 
For simplicity, let us assume that ESPs are distributed evenly over field space with some average 
inter-ESP distance y. In trapped inflation, the sensitivity of the inflationary trajectory to individual 
ESP encounters is greatly reduced: a change in ip* ± lowers the distance to about one half of the 
relevant ESPs, while the distance to the other half is raised. Since the contributions to the number 
of efoldings N from each ESP encounter are additive, the dependence of N on ip* ± vanishes in the 
limit of y — > 0. Thus, the effects computed in the previous sections are diminished for y smaller 
than some critical inter-ESP separation y cr it- 

Let us estimate this critical inter-ESP distance for a field-space of arbitrary dimensionality D. 
We assume the same quadratic potential for all inflatons, so that it has spherical symmetry. In the 
absence of ESPs, any trajectory from \(p\ = ip* ~ 15 to \(p\ = (fend ~ 1 yields the same inflationary 
dynamics. Thus, the total volume in which a slow-roll trajectory could lie is Vtot = {<f® — f^ndl^D 
where Vp = ir D / 2 /T(D/2 + 1) is the volume of a unit sphere in D dimensions. On the other hand, 
a grazing ESP encounter would only lead to observationally relevant contributions to observables if 
the ESP is within a distance of \i c ~ 10 -3 to a given inflationary trajectory. If ESPs are randomly 
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spread over field space, we can compute the probability that a single ESP in Vtot is within reach 
of the trajectory to 

PESP VbfcP-rfJ ' ' ' 

Thus, if N ESP ESPs are spread throughout Vtot, the probability that none of them is closely grazed 
becomes (1 — Pesp)^ esp - Demanding that the latter is less than 1/2 leads to the minimal number 
of ESPs in Vtot above which ESP encounters become likely, 

N™ = In (p ESP - . (80) 

The corresponding average inter-ESP distance is 



VcrU ~ V JV™? ^ J ■ ( 81 ) 

Note that 

Vcrit — >• in the limit of D — > oo. For D — 2, 10, 100 the minimal number 
of ESPs is N™™ ~ 1.7 x 10 4 , 10 38 , 10 413 with corresponding average inter-ESP distances of 
y cr it ~ 0.20,0.0030,0.00045. Thus, for low D, ESPs do not need to be particularly dense in 
order for a grazing ESP encounter to be likely. 

However, with growing dimensionality D, it becomes increasingly unlikely that only a few ESP 
encounters take place: either y > y cr n and no ESP is encountered, or y is (slightly) below y CT it and 



many ESPs are grazed. Thus, in the large D 



or trapped inflation at a terminal velocity 



26] 



imit we either expect undisturbed slow-roll inflaton, 



If trapped inflation takes place, the dominant correction, next to the ones stemming from a 
reduced speed, originate from backscattering of x-particles onto the inflaton condensate, which 



superpose to yield yet another near 
panied by non-Gaussianities j^. [40I- - 
been investigated in 



y 



scale invariant contribution to the power-spectrum, accom- 



421 ] . This superposition of bumps in the power-spectrum has 
411 ] for trapped inflation in one dimension. The corresponding computation in 
ligher dimensional field spaces for the power-spectrum and the three-point function is in progress 
60J 15 , but goes beyond the scope of the present article. 



To summarize, for inter-ESP distances of order y cr u and low D, several ESP encounters are 
likely, without leading to trapped inflation, and the results of this article are valid, namely, the 
dominant contribution to the power-spectrum can originate from individual ESP-encounters, nat- 
urally leading to observably large non-Gaussianities. 



15 Preliminary results show that the additional contribution to the power-spectrum acquires a blue tilt; hence, it 
cannot be the dominant contribution, but non-Gaussianities might still be observable. 
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IV. OBSERVATIONAL SIGNATURES FROM MODULATED TRAPPING 



We now consider the model of modulated trapping introduced in [8j. In this scenario, there is 
only one inflaton field, so that the inflaton trajectory is single-dimensional, but there exists another 
field, the modulaton a, which affects the strength of the trapping effect. As we will see below, the 
role of the modulaton on the primordial curvature perturbation is quite analogous to that of the 
perpendicular inflaton perturbation 5ip± discussed in the previous sections. 

More explicitly, the interaction between the inflaton (p and a bosonic field x is °f the form 

C = ~{m{a)-\{a) l p) 2 X \ (82) 
leading to a trapping event at time t p , with corresponding field value 

^ = W)- (83) 

Both the time and the strength of the trapping event depend on a, the modulaton field, which we 
assume to be light compared to the Hubble parameter during inflation. 

However, in contrast to we will not neglect its mass, and so its energy density after the 
trapping event should be taken into account. We will see that the modulatons potential in general 
also affects the spectral index and the scale dependence of /nl- If the modulaton field decays after 
the end of inflation, it acts as a curvaton field. In this section we assume that the energy density 
due to the modulaton during inflation is negligible (in order to be able to neglect its effect after 
inflation), so that 

£%, a) = V(ip) + W(a) ~ V(<p) ~ 3H 2 , V v ^>W, a e = e v ~ e H » e a . (84) 

However, the above relations leave the relative values of the second derivatives undetermined, so 
it is possible to have 1 3> \r) a \ > {r]^. 

We can write the primordial curvature perturbation as 



C(k) = £(k) + C G (k) + lUk) (C * C G ) (k) + ^g a {k) (C G * C * C) (k) , (85) 
where C G and C G are Gaussian, and it is a good approximation to treat the contribution from the 



fields as uncorrelated 



571 ] . The quantity f a , which is /nl in the limit that we neglect the inflaton 



field's perturbations (and similarly for <7tvl), is given by 



6 AN a<7 36 _ AN, aa(7 

5 U (AiY CT ) 2 ' 25 9a (AiV iCT ) 3 ' 1 j 
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where 



AN 



A 5 / 2 M 1 ' 2 



18vr 3 



H, 



^(2e p ) 1 /4 = AN{(Tp) 



(87) 



The power-spectrum is given by 

P f = P fv +Pc.=^(l-3)- 1 , (88) 
where we have defined S(fc) analogously to (fT7]) and the spectral index satisfies, (see ([IB])). 



n s - 1 = -(6 - 4H)e + 2(1 - H) % + 2Hr/, 
Then the non-Gaussianity and its scale dependence is given by 61 



62] 



"fNL = 



/nl 
dln|/ NL | 
5 In 



9NL 



4(1 -S)(2e + ryv + 
4(l-H)(2e + 77 CT -7^) + 



3 In A; 

-Mp , 



A AT, 



V 



(89) 



(90) 



(91) 



1 /6 



s U /7VL 

3V- 



(92) 
(93) 

. All terms, except those 



Note that S is an observable if we can observe the bi- and trispectrum 
involving AN, in the above equations should be evaluated at horizon crossing. rif NL is the sum of 
two independent parts: the first one is due to the presence of linear contributions from two scalar 
fields, as was discussed in Sec 2 of 62j, see in particular the mixed inflaton-curvaton scenario; this 



part vanishes if one neglects the perturbations due to the inflaton field (S = 1). The second part is 
due to non-linear interactions of the modulon field's perturbations, which vanish if the modulaton 
field has a quadratic potential, since the modulaton field perturbations obey a linear equation of 
motion (see Sec. 3 of [62|, in particular the interacting curvaton scenario). Observational prospects 



for nf NI are considered in 



63|, 



64]. 



In Sec. HH we considered two fields with equal mass and a quadratic potential (implying that 
all three first order slow-roll parameters are equal), so the scale dependence of /jvx is especially 
simple in this case, satisfying two independent consistency relations with the tensor-to-scalar ratio 
and spectral index 

r = 2n fNL = -4(n s " 1) = 16(1 - S)e* . (94) 

In the limit that perturbations come from the (initially) isocurvature field and H — > 1 all observ- 
ables, except for the three non-linearity parameters and the amplitude of the power-spectrum, are 
zero. 
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A. Special cases 



In [8| a special case was considered for which W(a) and all of its derivatives are negligible, with 
the consequence that 

n s -l = -(6-43)e + 2(l-5)7fc, (95) 
n fNL = 4(l-3)(2e-7fc). (96) 

If modulated trapping provides the dominant contribution to the power-spectrum, /nl is scale 
independent. Alternatively, if H is not too close to unity and e <C |t7 w |, a non-trivial consistency 
relation holds: 

n fNL = -2(n a - 1) . (97) 
In the special case of a modulaton dependent coupling, that is 

A = Id ' m = ~ 9G ' ( ' 98 ' ) 

9? p is independent of a. Then AiV oc o^ 2 and it follows that AN >a /AN >(Ta = 2a p /3. Assuming 
W{a) oc a q we find 

n fNL = = 4(1 " S)(2e + r? CT - ^) + |(g - 2)^ . (99) 

B. The modulaton as a curvaton 

After the trapping event the modulatons energy density needs to be taken into account. If 
a decays after the inflaton does, it can act as a curvaton, adding a third term to the curvature 
perturbation 

C = Cinf + Ctrap + Ccurv (100) 

Since Ctmp Ccurv oc for* the two latter terms are correlated, while Cinf oc 5ip* is uncorrelated 
to them. Unlike Ccurv ; Ctrap usually does not lead to any isocurvature perturbations since any 
particles created during the ESP encounter, and their decay products will be quickly diluted by 
the subsequent inflation, unless the encounter is close to the end of inflation. Due to the correlation 
properties of the three terms, we extend the definition of H to 
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where 



P c = P t ,_ t + P tt _ + P t 



Cinf 



Ct. 



Cc 



(102) 



Then all of the formula's given in (|89p - (|93|) for the spectral index and the non-linearity parameters 
remain valid. 

We would like to comment on a related idea: in the case of modulated reheating, the field which 
modulates the efficiency of reheating can also act as a curvaton [58|]. One could explicitly calculate 



both the adiabatic and isocurvature perturbations arisin g fr om this scenario, at the level of the 



power- and bispectrum using the formalism developed in 
using the formalism in 



651 ] , and at the level of the trispectrum 



6q which includes perturbations up to third order. We leave this project 



to future research. 



V. RELATING A GRAZING ESP ENCOUNTER AND MODULATED TRAPPING 

We have investigated two possibilities of how an ESP encounter can modulate perturbations 
in the presence of a second field: firstly, by putting the ESP in a two dimensional field space so 
that small fluctuations in the field perpendicular to the trajectory modulate the impact parameter 
and thus particle production, Sec. HH Secondly, by traversing the ESP head on, but including a 
dependence of the coupling to x and/or its bare mass on a second scalar field, Sec. IIVI How are 
these two scenarios related? 

The modulated trapping model is more general in the sense that the strength of the coupling and 
the position of the ESP is allowed to depend on the isocurvature field. However, it is significantly 
more restricted in the sense that a straight trajectory is assumed throughout inflation that traverses 
the ESP, i.e. the a field is not treated as an inflationary direction; consequently, no exponential 
suppression due to a non-zero impact parameter is possible (see e.g. the last term in ([35])). Hence 
neither scenario is a subset of the other, but they are similar in spirit and in certain cases may 
become effectively equivalent. 

Let us consider a modulaton-dependent coupling as in (|98p (the only explicit model considered 
in which can generate S — 1 with only one species of x particles), but also allow the model to be 
two-dimensional, that is we identify the perpendicular direction in Sec. [H] with the modulaton a. 
Then the change in the number of efoldings due to decreasing speed in field space in (j35[) becomes 



ANsl ~ % 4 e -^/(Mv^) ~ ^ W -2 e -3^/(M^i) (103) 

M 5 / 2 ^rri(p p M 5 /2 v ' 
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where we have chosen typical values ip p ~ 10 and m ~ 10 6 after the second ~ sign in order to 
simplify the analysis. This leads to 

5 a 3 ' 2 a 9 / 2 

We see that the contribution from the exponential term dominates unless a 3 < 10 _3 M, which is 
about the same condition for being able to neglect the exponential suppression in (| 103j) . Then in 
order to satisfy the COBE normalisation Vf ~ 10~ 10 (see (|65|) ). we need to satisfy a 3 H 2 10~ 4 /M 5 = 
10~ 10 ; hence, for this model to work in the regime of [81 we require M < 10 -2 and |c| < 10~ 2 . 
These conditions are consistent with those given in eq. (43) of [8J], but here we have also checked 
whether or not one can neglect the exponential suppression term. We see that there exists a 
consistent regime, but within a fine tuned range of initial values for a. 



A. Modulated trapping in a 2-dimensional field space 

The effect of the exponential suppression and fine tuning of the initial conditions is potentially 
reduced if we consider an Extra Species Locus (ESL) instead of an ESP. Two simple choices come 
to mind: a straight line ESL along ip = {p p , so that every trajectory passes through the ESP during 
inflation, or a compact ESL i.e. a circular one centred at the origin. Both choices are feasibly in 
moduli spaces, since they simply indicate that the ESL is independent of a field (the first choice) 
or an angular direction. We choose to focus on the second case only, without any theoretical bias 
either way. To make the two-dimensional field space explicit, we replace a by from here on. 
Thus, we wish to investigate an interaction Lagrangian of the form 

C = - 1 -g 2 (\^\-\^\) 2 X 2 , (105) 

with a modulated coupling, g 2 = <p 2 /M 2 . Then all trajectories become nearly straight lines and the 
exponential suppression term can be neglected for any trajectory. However, a simple calculation 
shows that approximately the same fine-tuning on the initial condition for <f2 is needed in order to 
have the correct amplitude of the power-spectrum if these perturbations dominate over the usual 
slow-roll contribution, i.e. IV2I ^ 10~ 2 is still required. 

If we had considered a constant coupling g, then all trajectories would be equally delayed by 
the ESL encounter and there would be no effect on the perturbation spectrum (in the case of a 
rotationally symmetric potential, which we are considering here). However, in general even an 
ESL with a constant coupling will generate perturbations provided that it breaks some symmetry 
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between trajectories with an isocurvature perturbation between them. For example, an elliptical 
ESL centred at the origin, or a circular ESL not centred at the origin, will both generate a con- 
tribution to the curvature perturbation at the time of particle production. Such a model would 
be similar to an inhomogeneous end of inflation 



67 



721 ]. which typically requires a very strongly 



7J- 



elliptical surface on which inflation ends as well as a fine tuning of initial conditions 

Why do we require this tuning, even when we are justified in neglecting the exponential sup- 
pression in (|103p ? A discussion of this point is given in the conclusion of [3]. In order for the 
perturbations associated with the trapping (</?2 direction) to dominate over those generated by the 
inflaton in the usual way (ipi direction), we require that 

/dN\ 2 ^/dN\ 2 1 > iq2 
\dip 2 J \d<pi) e* ~ 

The size of the isocurvature perturbations in each model is given by 5<p2 = 0(H) = 0(mi), 

independently of the background trajectory. Thus, these small perturbations correspond to a 

significant perturbation in the value of ip2 only if the background trajectory lies close to (f2 = 0, 

i.e. 5(f2/(f2 is non- negligible only for small tp2- This requirement is needed to satisfy the inequality 

(|106p for all of the models that we have considered so far. However, if <p 2 is too small, i.e. 

(f2 = 0(5(f2), then non-Gaussianities become large and the model is ruled out by observations 

16 



VI. CONCLUSIONS 



The possibility that the effective mass of a field coupled to the inflaton becomes small, leading 
to particle production, has been well studied in the case that this occurs either while modes of the 
scale observed today cross the horizon (a window of 5-10 efoldings) or shortly after inflation ends, 
i.e. preheating. However, there is a long period in between these two times (circa 50 efoldings) 
during which particle production may occur, and this case has been much less studied. Although 
we cannot directly observe any bump/feature this may cause on the power-spectrum on these 
scales, simply because we are unable to observationally probe these scales, we may see an effect on 
the usual CMB scales. 

In the simplest case, with only one inflaton field, there is no affect onto observables caused by 
such particle production, but if the strength or position of the ESP depends on a second field, 

16 In Sec. Mil the axis are chosen such that the ESP lies on the ipi-axis. The requirement that the value 992 is 
small is then identical to the requirement that the impact parameter jj, to the ESP is small in Sec. Mil Further, 
ip2 > 0(8ifi2) was written as ipj_ > 8ip QM which we imposed throughout Sec. HT1 and Sec. Mil 
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then this conclusion does not hold, as was studied by Langlois and Sorbo |8(. Here we focussed 
instead on multiple-field inflation. For most initial conditions, the trajectory in field space will 
not go through the ESP, but graze it, still leading to particle production and backreaction in the 
vicinity of the ESP. The dominant effect is usually the temporary slow down of the inflaton, as 
part of its kinetic energy is diverted to generating particles, which slow down the fields further via 
their backreaction; in addition the trajectory bends towards the ESP due to the attraction of the 
generated particles, whose mass grow as the inflaton moves further away. The expansion of the 
universe and associated dilution of these particles ensures that the effect only lasts a few efoldings 
before a slow-roll attractor trajectory is again reached. 

In multi-field scenarios the minimum distance between trajectories and the ESP varies due to 
the presence of isocurvature perturbations (those perpendicular to the velocity along the inflaton 
trajectory). The closer to the ESP that a trajectory comes, the greater the particle production 
and the larger the consequent slow-down of the trajectory hence the isocurvature perturbation 
is partially converted to the adiabatic perturbation on all scales which have already crossed the 
horizon. In order for this encounter to significantly affect the trajectory, it must pass close enough 
to be reasonably strongly distracted by the "beauty" of the ESP. For the models we have studied, 
an impact parameter of order 10 -3 Mp is required, which would represent a fine tuning of the 
initial conditions if only one ESP exists, but for an average ESP spacing of order 0.3Mp at least 
one such encounter becomes likely. An order of unity coupling constant between the inflaton and 
the produced particles is the most interesting case, so no tuning of this parameter is required. The 
case of a very dense set of ESP's, leading to multiple trapping events per efolding, was previously 
studied by D. and T. Battefeld in 261 ] . 



What are the most interesting consequences of such a grazing ESP encounter? Since this effect 
is asymmetric, most strongly affecting trajectories which pass closest to the ESP and most weakly 
those furthest away, this will generate a skewness in the probability distribution of trajectories, 
which is parametrised by Jnl (the amplitude of the bispectrum). The trispectrum is also affected, 
and the ratio t / {Q/nl > 1 provides a measure on which proportion of the observed linear 
perturbation where generated by the ESP encounter, with the inequality being saturated in the case 
that the encounter gives the dominant effect and the usual slow-roll perturbations are negligible. 
Since we have studied a quadratic potential with equal masses the scale dependence of the power 
spectrum and non-linearity parameters are suppressed, but for a more complex potential they 
would typically be a larger. 

We have also considered the possibility of a one dimensional enhanced symmetry locus, which 
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is a natural generalisation of the modulated trapping model to multiple fields. In this case every 
trajectory may pass through a point where the coupled field becomes massless, removing the need 
to fine-tune initial conditions, but one needs a strong difference in the strength of the particle 
production event between trajectories in order generate an observably large signature. 

In all cases we find observably large non-Gaussianities on CMB scales given that the power- 
spectrum is dominated by the effects brought forth via particle production during the intermediate 
50-efoldings of inflation. Thus, if PLANCK does not observe primordial non-Gaussianities, the 
interaction of inflatons throughout inflation are constrained. If inflatons are identified with moduli 
fields in string theory, we can thus probe not only the potential along the inflationary trajectory, 
but also the presence of enhanced symmetries in moduli space in the vicinity of the path the fields 
take. 



Acknowledgments 

T.B. and D.B. are grateful for hospitality at the University of Bielefeld and the Institude de As- 
troParticule et Cosmologie (APC, Paris). CB is grateful for hospitality received from the University 
of Goettingen. We would like to thank N. Barnaby for discussions. 



Appendix A: Conditions for instantaneous particle production 

We assume particle production to be instantaneous through the paper, but when is this a good 
approximation? If the non-adiabaticity parameter is large [l| 

— >1, Al 

where 



(t) ={k? + <?(<?- fy) 2 ) 1 '' , (A2) 



explosive particle production similar to pre- heating can take place. Evaluating the above at \<f(t) — 
<Pp\ = fi and focussing on k <C we find the condition 

< = /Vax , (A3) 

where v p = \<p p \- The subscript p denotes the time of particle production when the distance to the 
ESP is minimal. For larger impact parameters, the approximation in App. C of Q could be used 
(see Appendix [Bj, which provides an iterative procedure under the assumption that lo/lo 2 <C 1. 
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However, as we shall see in Sec. IIIIl the cases where an ESP encounter may be observable, but 
without observationally ruled out non-Gaussianities, indeed require [i < /i max anyhow. Then, the 
production event can be treated as instantaneous, as in Sec. [TTJ One may further check that the 
time-scale of particle production is considerably less than the Hubble time, so that the expansion 
of the universe and the inflationary potential can be neglected during this interval. 

Appendix B: Backreaction 

The backreaction of particle production onto the inflaton is a delicate point. This question is 
addressed in Q] for a two-dimensional scalar field space, without including a potential, and ignoring 
the cosmological evolution. 

In the present work, we follow the approximation proposed in jjj, which consists of decomposing 
the impact of particle production into two successive phases. In the first phase, which is very short, 
adiabaticity is strongly violated and one can compute, using the methods of preheating, the number 
of particles produced, assuming that the trajectory remains unperturbed. In the second phase, by 
contrast, the number of particles is assumed to be conserved (i.e. simply diluted by the volume 
expansion in our cosmological case) and the backreaction on the trajectory is taken into account by 
considering the interaction between the scalar field and the produced particles, or more specifically 
the energy density of the particles, given by 



where n is the number density of particles and m x {(p) = g\(p — <p p \ their effective mass. 

For the model discussed in this paper, the dynamical evolution of the inflaton is thus described 
by the Lagrangian 



Vint = p x (<P) = g\ ( P-0p\n 



(Bl) 




(B2) 



with 




(B3) 



The associated equation of motion is given by 



" } _^ 2 — * 

f + 3Hif + m <f 




Op\ 3 <f_ ~ <PeSP 



(B4) 
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a J \\tp-<p ESP 
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If one decomposes this equation along the directions parallel and orthogonal to the initial 
(unperturbed) trajectory, one finds 

n - 



<p\\+3H<p\\ +mV|| = -gnJ^) 6(t-fr)-= J y||ESP = , (B5) 

^/ (^|| - ^||ESP) 2 + - 

<p ± + 3H<p ± +m 2 ip ± = -gnJ^-Y Q(t-t p ) ip± ~ f * , (B6) 

° yJ(<P\\ -^||esp) 2 + (V± -V) 2 

where we have used the fact that the gradient of the potential vanishes in the orthogonal direction 
and <^_i_esp = f 1 - (We put the ESP above the (^y-axis here, opposite to Fig. Q]and Fig. [2] but in line 
with ([48 P and the numerical plots in this paper). 

The numerical results presented in the main body of this paper are based on the above system of 
equations. Moreover, the approximate equations (j29p and ([48D . which are integrated analytically, 
are derived from ([B5P and (jB6j) by assuming (pj_ <C /ti. 

When the above approximation is not valid, in particular when particle production cannot 
be considered as instantaneous, one must in general solve numerically the system involving the 
evolution of the mode functions of the x field combined with the equation of motion for the scalar 
field, as done in ((]] in the case of a scalar field with V = in Minkowski spacetime 17 . 
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